Published by Institute of Physics Publishing for SISSA/ISAS 

Received: 
Revised: 
Accepted: 



Neutralino Dark Matter in Mirage Mediation 



Kiwoon Choi, Kang Young Lee, Yasuhiro Shimizu 

Department of Physics, KAIST, Daejeon 305-701, Korea 



E-mail: kchoiOmuon .kaist . ac .kr , kylee@muon.kaist.ac.kr 



shimizuOmuon . kaist . ac . kr 



Yeong Gyun Kim 

Astrophysical Research Center for the Structure and Evolution of the Cosmos, 
Sejong University, Seoul 14-3-747, Korea 



E-mail: ygkimOmuon. kaist . ac .kr 



Ken-ichi Okumura 

Department of Physics, Kyushu University, Fukuoka 812-8581, Japan 



E-mail: okumuraQhiggs . phys . kyushu-u .ac.jp 



Abstract: We study the phenomenology of neutrahno dark matter (DM) in mirage me- 
diation scenario of supersymmetry breaking which results from the moduli stabilization in 
some string/brane models. Depending upon the model parameters, especially the anomaly 
to modulus mediation ratio determined by the moduli stabilization mechanism, the na- 
ture of the lightest supersymmetric particle (LSP) changes from Bino-like neutralino to 
Higgsino-like one via Bino-Higgsino mixing region. For the Bino-like LSP, the standard 
thermal production mechanism can give a right amount of relic DM density through the 
stop/stau- neutralino coannihilation or the pseudo-scalar Higgs resonance process. We also 
examine the prospect of direct and indirect DM detection in various parameter regions of 
mirage mediation. Neutralino DM in galactic halo might be detected by near future direct 
detection experiments in the case of Bino-Higgsino mixed LSP. The gamma ray flux from 
Galactic Center might be detectable also if the DM density profile takes a cuspy shape. 
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1. Introduction 

Low energy supersymmetry (SUSY) is a promising candidate for physics beyond the stan- 
dard model at TeV scale In addition to solving the naturalness problem associated 
with the weak to Planck scale hierarchy M^^ak/Mpi ~ 10~^^, low energy SUSY provide a 
natural candidate for the cold dark matter (DM) in the universe the neutralino being 
the lightest supersymmetric particle (LSP), under the assumption of ii-parity conserva- 
tion. The assumption that the LSP neutralino constitutes DM constrains SUSY model in 
various ways. The first constraint comes from the requirement that the LSP of the model 
should be a neutralino. Although not severe, the current limits on the direct or indirect 
detection of DM also provide an additional constraint on the model. Finally, if the mech- 
anism of cosmological neutralino production is specified, the model parameters should be 
in the range giving the correct DM abundance measured by WMAP 1^]: 

0.085 < noAih^ < 0.119. (1.1) 

Recent development in string moduli stabilization [Q, || has led to a new pattern of 
soft SUSY breaking terms which has not been explored before [^. In KKLT-type moduli 
stabilization scenario |5|, the volume modulus T is stabilized at SUSY AdS vacuum by 
non-perturbative effect, and this SUSY AdS vacuum of T is lifted to a phenomenologically 
viable dS or Minkowski vacuum by SUSY breaking brane. Such set-up leads to a mass 
pattern @: 



where mx and denote the modulus mass and F-component, respectively. If the SUSY 
breaking brane is sequestered from the visible sector, the soft masses of visible fields are 
determined by the modulus mediation |^ of 0{F'^ /T) and the anomaly mediation Q of 
O (7713/2 /Svr^) which are comparable to each other if the gravitino mass 777-3/2 ~ 10 TeV 
as required to give the weak scale size of soft masses. An interesting consequence of this 
mixed mediation is that soft masses are unified at a mirage messenger scale Q 

M^i, = Mct/T ^ , (1.3) 

where a is a parameter of order unity which represents the anomaly to modulus mediation 



ratio^. This feature named as mirage mediation |10| is due to the particular property of 
anomaly mediation which is closely related to the renormalization group (RG) evolution 
of soft parameters. KKLT-type moduli stabilization typically gives a positive a = 0{1), 
and thus a mirage messenger scale hierarchically lower than Mqjjt- The corresponding 
low energy superparticle masses have a quite different pattern from those in other SUSY 



breaking scenarios such as mSUGRA, gauge mediation and anomaly mediation |9L h]l, 12 



i3i. 



In this paper, we wish to study some phenomenological aspects of neutralino DM in 
mirage mediation, particularly examine the parameter range giving LSP neutralino, the 
relic DM density under the assumption of standard thermal production, and the prospect 
of direct or indirect detection of DM neutralino. As we will see, when the anomaly to 
modulus mediation ratio a increases from zero to a positive value of order unity, or equiv- 
alently the mirage messenger scale Mmir = Mgc/^ (7773/2 /Mp;)"/^ decreases from Mqut 
to a hierarchically lower scale, the nature of LSP neutralino is changed from Bino-like to 
Higgsino-like via Bino-Higgsino mixing region. The enhanced Higgsino component can be 
understood by the gluino, Wino and Bino mass ratios in mirage mediation: 

M3 : M2 : Ml ~ (1 - Q2a)gl : (1 + 0.1a)gl : (1 + OMa)gl, (1.4) 

where a ^ represents the effects of anomaly mediation. Thus, compared to the mSUGRA- 
type pure modulus mediation (a = 0), the gluino to Bino mass ratio for a ~ 1 is signif- 
icantly reduced. This results in smaller Irnj^^l/Mf at the weak scale, and thus a smaller 
Higgsino to Bino mass ratio when the electroweak symmetry breaking condition is imposed. 
(Here rnj^^ is the soft SUSY breaking mass-square of the up-type Higgs doublet Hu-) As 
a consequence, for q ~ 1, the LSP neutralino has a sizable Higgsino component over a 
large fraction of the parameter space, and this Higgsino component eventually becomes 
dominant for a larger value of a. For the parameter region leading to Bino-like LSP, the 
thermal production mechanism can give a right amount of relic DM density through the 
stop/stau- neutralino coannihilation process or the pseudo-scalar Higgs resonance effect. 



^Here we are following the notations of Ref. [^. We warn the readers that some subsequent works [|l^, 
are using a differently defined a which corresponds to a^^^ — 167r'^/[aours ln(Afp;/m3/2)l — 4.9/aours- 
We stress that our definition of a is more convenient for rnatching the mirage mediation to underlying 
moduli stabilization model. For instance, the original KKLT moduli stabilization gives Ofours ~ 1, and 
many of its generalizations give a rational value of Qours- 



depending upon the value of tan/3 = (Hu) / {H^) . In overall, compared to mSUGRA sce- 
nario, a significantly larger fraction of the parameter space of mirage mediation can give 
the WMAP DM density under the assumption of the conventional thermal production of 
neutralino LSP, while satisfying all known phenomenological constraints. As for the DM 
detection, neutralino DM in galactic halo might be detected by near future direct detec- 
tion experiments in the case of Bino-Higgsino mixed LSP which can be obtained over a 
large fraction of the parameter space. The gamma ray flux from galactic center might be 
detectable also if the DM density profile takes a cuspy shape. 

Some phenomenological and cosmological aspects of mirage mediation have been inves- 



tigated by several groups |ri| |13|, |15|. In particular, the properties of neutralino 



DM have been studied in Ref. The sub-GUT CMSSM model of Ref. |T6| also shares 
a common feature with mirage mediation as it assumes the unification scale of soft pa- 
rameters below the GUT scale. In this paper, we provide a more extensive analysis of 
neutralino DM for certain range of model parameters which are expected to be obtained 
in KKLT-type moduli stabilization scenarios together with a more detailed study of the 
direct and indirect DM detections. Our results agree qualitatively well with Ref. |13| when 
the considered parameter space overlaps. 

Recent study of moduli cosmology in KKLT-type moduli stabilization indicates 
that if the early universe underwent a period dominated by the coherent oscillation of 
moduli, there can be too many gravitinos and neutralino LSPs produced by moduli decays, 
which would spoil the Big-Bang nucleosynthesis or overdose the universe. This requires a 
mechanism to dilute the primordial moduli oscillation, e.g. a thermal inflation as proposed 
in |17] . In the presence of such thermal inflation, the neutralino DM can be produced either 
by the conventional thermal production mechanism or by the decays of flaton triggering 
thermal inflation, depending upon the details of reheat procedure. Throughout this paper, 
we will focus on the relic density of neutralino DM produced by the conventional thermal 
production mechanism while leaving an additional non-thermal production ^] as an 
open possibility. 

In fact, if the anomaly to modulus mediation ratio a takes a value giving a mirage 
messenger scale Mmir = Mcf/j^ (171,3/2 /Afp;)°/^ much lower than the intermediate scale, the 
squark/slepton mass-squares have negative values at certain /iz^/i renormalization point > 



/ic, although they become positive at lower renormalization point around TeV [^, 12 1. As 
long as the low energy squark/slepton mass-squares are positive'^, the model has a correct 
color /charge preserving (but electroweak symmetry breaking) vacuum with squark/slepton 
fields (j) = 0. Still, tachyonic squark/slepton mass-squares at high renormalization point ^ > 
fic might give a deeper color /charge breaking (CCB) minimum or a unbounded from below 
(UFB) direction in the effective potential at (p > fic pOC , which would make the color /charge 
conserving vacuum metastable. In the presence of such CCB minimum or UFB direction, 
there are two points to be clarified to make sure that the model is phenomenologically 
viable. One first needs a cosmological scenario which allows our universe to be settled down 
at the correct vacuum, not at the CCB minimum or UFB direction. The second is that the 



^It has been argued that a negative stop mass-square at high renormalization point is helpful for ame- 
liorating the fine-tuning for the electroweak symmetry breaking in the MSSM llSf. 



color/charge preserving vacuum should be stable enough against the tunnelling into CCB 
minimum or UFB direction. The first point might depend on the detailed history of the 
early universe. However in view of that squark/sleptons get large positive mass-squares 
in the high temperature limit, it is a rather plausible assumption that squark/sleptons 
are settled down at the color /charge preserving minimum after the inflation [^]. As for 
the vacuum stability, it has been noticed that the potential barrier between (p = and 
(j) ^ fj-c gives a tunnelling rate much less than the Hubble expansion rate as long as ^ 10 
TeV |22]. In mirage mediation with positive a = 0{1), once one requires a successful 
electroweak symmetry breaking, is higher than 10^ TeV, thus satisfies safely the stability 
condition. In this paper, we do not take the existence of CCB minimum or UFB direction 
at large squark/slepton value \(/)\ ^ 10 TeV as a serious problem of the model as long as 
a good color /charge preserving and electroweak symmetry breaking vacuum exists, and 
focus on the phenomenology of the model under the assumption that we are living in the 
color /charge preserving local minimum which is stable enough to have a lifetime longer 
than the age of the universe. 

The organization of this paper is as follows. In section II, we review the basic features of 
mirage mediation. In section III, we examine the prospect of neutralino DM in intermediate 
scale mirage mediation scenario for several different choices of the matter/Higgs modular 
weights. In section IV, we extend the analysis to general values of the mirage messenger 
scale. Section V is the conclusion, and Appendix A contains a summary of our convention 
and notation. 



2. Mirage mediation 



Mirage mediation is a natural consequence of the KKLT-type moduli stabilization scenario 
satisfying the following two assumptions: (i) the modulus T (or dilaton) which determines 
the standard model gauge couplings is stabilized by non-perturbative effects and (ii) SUSY 
is broken by a brane-localized source which is sequestered from the visible sector. A well 
known example of such set-up is the KKLT moduli stabilization H] in type IIB string 
theory^ A similar but simpler example would be 5D brane model with a flat interval in one 
side and an warped interval in other side, in which SUSY breaking brane is introduced at 
the IR fixed point of the warped interval. 

Under these two assumptions, the 4D effective action of the visible sector fields and 
the gauge coupling modulus T is given by 



-3CC*e 



-K/3 



2a2 



+ h.c. 



(2.1) 



where C = Cq + F^O"^ is the chiral compensator superfield, fa are the holomorphic gauge ki- 
netic function of the visible sector gauge fields, K and W are the effective Kahler potential 



^As was noticed in due to the effect of throat vector multiplet, the sequestering might not be 

precise enough in the case of KKLT compactification of type IIB string theory. The size of non-sequestered 
soft scalar mass induced by the exchange of throat vector multiplet is quite sensitive to the unknown details 
of compactification, however there exists a reasonable parameter limit in which the non-sequestered effects 
can be safely ignored [E3|. 



and superpotential of the visible matter superfields <I>j and the gauge coupling modulus T, 
which would be obtained by integrating out heavy moduli. As long as the SUSY breaking 
brane is sequestered from the visible gauge and matter superfields, its low energy conse- 
quence can be described by a simple spurion operator of the form T'mtO'^O'^, independently of 
the detailed dynamics on the SUSY-breaking brane. Assuming an axionic shift symmetry: 

C/(1)t : Im(r) Im(r) + real constant (2.2) 

which is broken by non-perturbative term in the superpotential, the model is given by 

K = Ko{T + T*) + Zi{T + r*)^>*$i. 



W = w- Ae-"^ + ^Xijk<^i<^j<^k, 
fa = kT + A/„ 

rmt = Vm{T + T*), (2.3) 



where a and k are (discrete) real constants, while Afa,w,A and Xijk are complex effective 
constants obtained after heavy moduli are integrated out. The axionic symmetry ( |2.2| ) 
ensures that Kq, Zi and Vmi depend only on T + T*, and a and k are real constant. With 
these features, the resulting gaugino masses and trilinear A parameters preserve CP as was 



pointed out in [125 1. 



There might be various ways to generate the modulus superpotential Wo = w — Ae~°''^ 
stabilizing T. Generically the non-perturbative term e~"^ can be induced by either a 
gaugino condensation of T-dependent hidden gauge interaction or a stringy instanton whose 
Euclidean action is controlled by T. As for the constant term w, it might be induced by a 
fine-tuned configuration of fluxes as in the original KKLT scenario ||5|, or alternatively by 
T-independent non-perturbative effect whose strength is controlled by heavy moduli p6| . 
As we will see, the non-perturbative stabilization of T by Wq generates a little hierarchy 
between the modulus mass and the gravitino mass: 

^ ~ aT ~ ln{Mpi/my2), (2.4) 
which results in a little suppression of the modulus F-component: 

^3/2 ^3/2 



ml/2 ^3 
T niT \n{M PI / 7713/2)' 



(2.5) 



Then, F'^ /T naturally has a size comparable to the anomaly mediated soft mass of 0(7713/2 /4'7r^ 
for the gravitino mass around TeV. If the SUSY breaking source is sequestered from the 
visible sector, the soft terms of visible fields are determined by the modulus mediation of 
0{F'^ /T) and the anomaly mediation of O (7773/2 /4'7r-^) which are comparable to each other. 
This leads to a mirage unification Q of soft masses at a scale hierarchically different from 
the gauge coupling unification scale Mqut- 

In the Einstein frame, the modulus potential from (|2.1|) is given by 



Wot = [{dTdfKo)-^\DTWo\^ - 3|t^o|'] + Mift, (2.6) 



where 



Wo = w- Ae-'^^, Mift = e2^»/3piift. (2.7) 



The superspace lagrangian density (2.1) also determines the auxiliary components of C 
and T as 

pT = _e-^o/2 (a^a^.i^o)-! {DtW^Y , (2.8) 

where = c^^^'^Wq. Note that one can always make both w and A real by appropriate 
U{l)ji and U{1)t transformations. In such field basis, the U{1)t invariance of Kq assures 
that both ^^'^ real. In the following, we will use this field basis in which the 

CP invariance of soft parameters is easier to be recognized. 

To stabilize T at a reasonably large value while having 777-3/2 hierarchically smaller 
than Mpi, one needs to assume that w is hierarchically smaller than A in the unit with 
Mpi = 1. Since w ~ 7773/2 and one needs 7)73/2 ~ 10 TeV to get the weak scale superparticle 
masses, ln{A/w) typically has a value of 0(471^). It is then straightforward to compute 
the vacuum values of T and by minimizing the modulus potential (|2.6| ) under the fine 
tuning condition (Vtot) = 0. At leading order in e = 1/ ln(j4/7u), one finds § 

aT ~ ln{A/w), 

^ L 39t In(Piift) 



T + T* ln(^/7i;) V 29tKo 

which shows that F'^ /T is indeed of the order of 7713/2/471^ for \\i{A/w) ~ 47r^. 

Let us consider the soft terms of canonically normalized visible fields derived from the 
4D effective action (^.ll): 



Aoft = -^MaA"A" - ]^mf\ct)i\^ - ^Aijkyijk(t>i(t>j(l}k + h.c, (2.10) 

where A*^ are gauginos, (pi are the scalar component of $j and yijk are the canonically 
normalized Yukawa couplings: 

y^^ = ^^^-fe (2.11) 

For F^^ jT ~ 7n3/2/47r^, the soft parameters at energy scale just below Mqut are deter- 
mined by the modulus-mediated and anomaly-mediated contributions which are compara- 
ble to each other. One then finds [ptl 



Ma = Mo + -^bagl 



167r2 

"^3/2 



Aijk = Aijk - (7i + 7j + 7fc), 

,2 _ ~ 2 ^m.,.^2 



where Mq, Aij^ and fhi arc the pure modulus-mediated gaugino mass, trihnear ^-parameters 
and sfermion masses which are given by 

Mo = F'^9rlnRe(/„), 

= -F^F^*dTdf\n{e-'''^'^Zi), 



e-KoZ,ZjZk 



= Ai + Aj + Ak for Ai = F^Ot ln{e-^°/^Zi). 



(2.13) 



Here we have used that the holomorphic Yukawa coupHngs A^-fc arc T-indcpcndcnt constants 
as a consequence of the axionic shift symmetry U{1)t, and the one-loop beta function 
coefficient ba, the anomalous dimension 7^ and its derivative 7^, and 9i are defined as 



ba = -3tr (r2(Adj)) +J2iT (r2(^,)) , 

i 

7i = 2^5^C2»(</.,)-^^|y,X 

a 

2 dli 



jk 



1; 



Sir' 



c/ln/Lt' 



jk 



2 -^ijk 

Mo"' 



(2.14) 



where the quadratic Casimir C|(0j) = (A^^ — 1)/2N for a fundamental representation (f)i of 
the gauge group SU{N), C^{(pi) = qf for the U{1) charge qi of (f)i, and ujij = YlklVikW^kl 
is assumed to be diagonal. In Appendix A, we provide a summary of our convention and 
notation. 

For our later discussion, it is convenient to define 



m 



a = 



i/2 



Mo ln(Mp,/m3/2) 



a,: — 



A 

Mo' 



C,; 



- 2 

Mf' 



(2.15) 



where a represents the anomaly to modulus mediation ratio, while ai and Cj parameterize 
the pattern of the pure modulus mediated soft masses. Then the boundary values of soft 
parameters at Mqut are given by 



A 



ijk 
2 



Mo 

Mr 



m: 



M, 



+ ^-9a<^ 

(ai + ttj + ak) j^g^2 + + 



ln(Mp;/m3/2) 
^ 16^^ 



ln(Mp;/m3/2 
167r2 



where 



= 45Z5»^2 (<^i) - IVijkWai + aj + ak). 

a jk 



(2.16) 



(2.17) 



In this prescription, generic mirage mediation is parameterized by 

Mo, a, Ui, Ci, tan/3, 



(2.18) 



where we have replaced the Higgs mass parameters /i and B by tan /? = {Hu) / {H^) and 
Mz as usual. As we will see, this parameterization of mirage mediation is particularly 
convenient when one compute the mirage mediation parameters from underlying SUGRA 
model. For instance, a, ai and q are given by simple rational numbers in minimal KKLT- 
type moduli stabilization. 

Taking into account the 1-loop RG evolution, the soft masses of (l2l^ ) at Mgut leads 
to low energy soft masses which can be described in terms of the mirage messenger scale: 

^^^^ (2.19) 



For instance, the low energy gaugino masses are given by Q 



Ma{yL) = Mo 



1 . 2 / \ 1 I -^mir 



rMa(M) In 



:Mo, 



(2.20) 



showing that the gaugino masses are unified at Mmir, while the gauge couplings are unified 
at Mgut- The low energy values of Aijk and mf generically depend on the associated 
Yukawa couplings i/ijk- However if Uiji. are small enough or 



ak = Ci + Cj + Cfc = 1 for yijk ~ 1, 



(2.21) 



their low energy values are given by 
Mo 



Mi 



ai + aj + Ofc + ^iliifJ') + IjifJ-) + 7fe(/w)) In 



Mrr 



If 1 djiifl) ( Afmir 



In 



(2.22) 



where Yi is the f7(l)y charge of (pi. Quite often, the modulus-mediated squark and slepton 
masses have a common value, i.e. Cq = c^. Then, according to the above expression of low 
energy sfermion mass, the 1st and 2nd generation squark and slepton masses are unified 
again at the mirage messenger scale Mmir- 

In KKLT compactification of type IIB string theory [|[, T corresponds to the Calabi- 
Yau volume modulus, and the uplifting brane is located at the end of warped throat. In 
this case, the uplifting operator is sequestered from T M: 



lift 



0. 



(2.23) 



Also, the minimal KKLT compactification of IIB theory gives 
Ko = -3ln{T + T*), Z, ^ 
fa = kT, Wo = w 



_)_ 

Ae~"^ {A = 0(1)), 



(2.24) 



o 
o 

o 
o 

o 
o 
o 



-8- 



where the modular weight rii is a rational number depending on the origin of matter 
super field <I>j. Then from (2.13) and (|2.9| ), one immediately finds 



a = 1, ai = Ci = 1 - Ui, (2.25) 

giving an intermediate mirage messenger scale: 

Mmir ~ 3 X 10^ GeV. (2.26) 

If <I>j lives on the entire world-volume of D7 brane from which the visible gauge bosons 
originate, the corresponding modular weight rii = 0. However, if $j is confined in the 
intersections of D7 branes, rii has a positive value, e.g. nj = 1/2 or 1. In Fig. ||, we depict 
the RG evolution of gauge couplings and soft masses in intermediate scale mirage mediation 
scenario with rii = 0, i.e. a = Oj = q = 1, which shows that indeed the gaugino masses and 
the 1st and 2nd generations of squarks and slepton masses are unified at Mmir ~ 3 x 10^ 
GeV as indicated by ( 2.20 ) and ( p. 22 ). Note that still the gauge couplings are unified at the 



conventional GUT scale Mqut 2 x 10^^ GeV. In view of its minimality, intermediate scale 
mirage mediation can be considered as a benchmark scenario, thus we perform a detailed 
analysis of neutralino DM in intermediate scale mirage mediation in the next section. 

It is in fact easy to generalize the compactification to get different values of the mirage 
mediation parameters a, a, and q. For instance, the compactification can be generalized 
to have a dilaton- modulus mixing in gauge kinetic functions ||2^, ^ and/or the non- 
perturbative superpotential |29|. For the case of type IIB compactification, a nonzero 



gauge flux on D7 branes can generate such dilaton-modulus mixing in D7 gauge kinetic 



functions |28], which would result in 

fa = kT + ISo, Wo = w- ^e-("^+''^o) {A = 0(1)), (2.27) 

where denotes the vacuum value of the string dilaton S which is assumed to get super- 
heavy mass from RR and NS-NS 3-form fluxes, and k,l,a and b are real parameters. In 
such IIB compactiflcations, the coeflicients of T in D7 gauge kinetic functions and non- 
perturbative superpotential, i.e. k and a, are positive, however the coefficients of S, i.e. / 
and b, might have both signs under the conditions: 

fcRe(T) + lRe{So) ~ ^ 



2 

9gut 



aRe(T) + 6Re(S'o) ~ \n{Mpi/m^i2) ~ 47r2. (2.28) 

Assuming that e~^°/^Zi has the same form as the minimal model (2.24), it is then straight- 
forward to find that the mirage mediation parameters are given by 

l + Ri 

(l+i?2)(l+i?3)' 
(l-n,)(l+i?i), 

{l-ni){l+Rif (2.29) 



a 



2.5 p , , , f , , , , , , , f ,— r 




Figure 1: RG evolution of (a) gauge couplings a^, (b) gaugino masses Ma, (c) sfermion and Higgs 
masses rm, (d) trilinear A parameters in intermediate scale mirage mediation with ai — Ci = I. 
Here we choose Mq = 800 GeV and tan/3 = 10. 



where 

Again, for Vuft induced by an uplifting brane at the end of warped throat, R3 = 0. However, 
1 + i?i and 1 + R2 can have a variety of (positive) values under the condition ( |2.28| ). As 
a result, the anomaly to modulus mediation ratio a can easily have any (positive) value 
within the range of order unity. 

A particularly interesting model is the TeV scale mirage mediation with 

a = 2, = = 0, agg + au^ = cq^ + = 1, (2.31) 

where Q3 and C/3 denote the left-handed and right handed top-quark superfields, respec- 
tively. This particular model has been claimed to minimize the fine tuning for the elec- 
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Figure 2: RG evolution of (a) gaugino masses and (b) sferniion and Higgs masses in TeV scale 
mirage mediation. Here we fixed Mq — 800 GeV, ajj = cjj = and am = cm — 1/2. 



troweak symmetry breaking in the MSSM |3C]. In Fig. g, we depict the RG evolution of soft 



masses in TeV scale mirage mediation model with an = ch = and om = cm = 1/2, where 
the subscripts H and M stands for the MSSM Higgs doublets Hu^d and the quark/lepton 
matter superfields, respectively. 

Note that the squark/slepton mass-squares renormahzed at high energy scale, e.g. 
at a scale near Mqut^ are negative in this model, while the values at low energy scale 
below 10^ GeV are positive. In fact, such tachyonic high energy squark/slepton mass- 
squares is a generic feature of mirage mediation for a > ac where the precise value of Qc 
depends on and q, but not significantly bigger than 1 in most cases. As long as the 
low energy squark/slepton mass-squares are positive, the model has a correct color/charge 
preserving (but electroweak symmetry breaking) vacuum. For instance, the TeV scale 
mirage mediation model of Fig. ^ has a such vacuum which is a local minimum of the 
scalar potential over the squark/slepton values |(/)| < 10^ GeV. On the other hand, tachyonic 
squark mass-squares at the RG point /i > 10^ GeV indicates that there might be a deeper 
CCB minimum color/charge breaking (CCB) or an unbounded from below (UFB) direction 
at > 10^ GeV. One then needs a cosmological scenario which allows our universe to be 
settled down at the correct vacuum with (/> = 0. In view of that the squarks and sleptons get 
large positive mass-squares in the high temperature limit, it is rather plausible assumption 
that squark/sleptons are settled down at the color /charge preserving minimum after the 
inflation One still needs to confirm that the color /charge preserving vacuum is stable 
enough against the decay into CCB vacuum. It has been noticed that the corresponding 
tunnelling rate is small enough, i.e. less than the Hubble expansion rate, as long as the 



RG points of vanishing squark/slepton mass-squares are all higher than lO'' GeV |2|, |l| 
which is satisfied safely by the TeV scale mirage mediation of Fig. 0. 



3. Neutralino DM in intermediate scale mirage mediation 



In this section, we examine the prospect of neutrahno DM in intermediate scale mirage 
mediation scenario. As was noticed in the previous section, the minimal KKLT-type model 



( 2.24 ) with a sequestered uplifting brane gives a = 1, thus an intermediate mirage messen- 



ger scale 

Mmir ~ MGUT{ms/2/Mpi)^/^ ~ 3 X 10^ GeV. (3.1) 

In this minimal set-up, the discrete parameters and describing the modulus mediated 
A-parameters and sfermion masses are determined to be = Cj = 1 — nj, where nj denote 
the modular weights of matter and Higgs superfields. Throughout this paper, we will 
assume = Cj and consider the following four different cases: 

(aH = CH, aM = CM) = (l, 1), (0,1), (0, i). (3.2) 

We also choose the Higgsino mass parameter ^ > in light of the experimental value of 
the muon anomalous magnetic moment which favors positive /i [^l|, and treat tan/3 = 
{Hu)/{Hd) as a free parameter without specifying the origin of the corresponding /x and B 
parameters. We then obtain the parameter range of the model for which the LSP is the 
lightest neutralino as well as the relic neutralino DM abundance under the assumption of 
thermal production, and finally the direct and indirect detection rates of the neutralino 



LSP using the DarkSUSY routine p2| 



3.1 Parameter region of neutralino LSP and thermal relic density 

In Fig. ^ we show the neutralino LSP region and the thermal neutralino relic density in 
intermediate scale mirage mediation scenario on the (tan/3, Mo)-plane for the values of Oj 
and Ci specified in Eq. ( |3.2| ). We computed the sparticle mass spectrum at the electroweak 
scale by solving the RG equations with the boundary condition (2.12) at Mqut- Our 



results show that in all cases there is a large parameter region for which the LSP is given 
by the lightest neutralino. 

Fig. ^.a is the result for the case in which Oj = q = 1 for both matter and Higgs 
multiplets. In this case, large tan/3 > 34 (grey color) for which the tau Yukawa couplings 
becomes sizable gives stau LSP (see also Fig. ^.a), while small Mq < 0.5 TeV (green color) 
gives stop LSP. In the remaining region, the LSP is the lightest neutralino which turns out 
to be Bino-like. Small tan/3 ^ 3 is excluded by the Higgs mass limit m/j > 114 GeV. Under 
the assumption that the DM neutralinos are produced purely by the conventional thermal 
production mechanism, the magenta stripe corresponds to the parameter region giving a 
relic neutralino density consistent with the recent WMAP observation |^]: 

0.085 < J)dm/i^ < 0.119 {2a level). (3.3) 

In the region below the magenta stripe, O^/i^ < 0.085, while ^^^^^ ^ 0.119 for the upper 
region. Thus the (cyan) region below the magenta stripe (but above the stop LSP region) 
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Figure 3: Parameter region of neutralino LSP and the thermal relic density depicted on the plane 
of (tan/3, il/o) in intermediate scale mirage mediation models with {ai,Ci) specified in Eq. (3.2). 



can be phenomenologically viable if additional DM neutralinos were produced by non- 
thermal mechanism such as the decays of flaton in thermal inflation |17]. 

Although the neutralino LSP is Bino-like in this particular intermediate scale mirage 
mediation, the WMAP mass density is obtained for a rather heavy neutralino mass m^o > 
450 GeV. This can be understood by Fig. ^a which shows the masses of the lightest 
neutralino, lighter stop and stau, and also the pseudo-scalar Higgs boson as a function of 
tan/3 for the model with q = 1, Oj = Cj = 1 and AIq = 800 GeV. Around tan/? ~ 20, 
the pseudoscalar Higgs mass becomes same as 2m^o, leading to a resonant enhancement 
of neutralino annihilation through the s-channel pseudo-scalar Higgs exchange. For other 
values of tan /?, the neutralino mass is somewhat close to the stop mass (or to the stau 
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Figure 4: Particle masses as a function of tan/3 in intermediate scale mirage mediation models 
with (flijCi) specified in Eq. (3.2). Here, we fixed Mo = 800 GeV. 



mass at tan/3 ~ 34), making the stop-neutralino or stau-neutralino coannihilation process 
becomes efficient. In Fig. ^, we depicted ^i^/i^ as a function of tan/3 for Mq = 800 GeV, 
which shows clearly the effect of Higgs resonance at tan /3 ~ 20 and also the effect of 
stop/stau coannihilation effects for other values of tan/3. On the plot, the dotted line 
corresponds to the relic density computed without including coannihilation effects. It 
indicates that the stop/stau-neutralino coannihilation plays a crucial role for the relic 
neutralino density to have the WMAP value (U) for Mq = 700 ~ 800 GeV and tan /3 
outside the Higgs resonance region. Note that in intermediate scale mirage mediation with 
flj = Cj = 1) the pseudoscalar Higgs resonance condition mA — 2?t),^ is satisfied for smaller 
value of tan/3 compared to the mSUGRA case. This can be understood by noting that 
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the low energy gaugino masses in mirage mediation are more compressed compared to 
mSUGRA, e.g. M3/M1 ~ 2.3 in the intermediate scale mirage mediation with a = 1, 
while M3/M1 ~ 6 in mSUGRA. For a given value of Mi, smaller M3 gives smaller /i and 
~ rn'jj^ + ^u^ at the electroweak scale, thus the pseudoscalar Higgs resonance appears 
at smaller value of tan/3 compared to mSUGRA case. 
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Figure 5: fl^h as a function of tan/? in intermediate scale mirage mediation with ai — Ci — 1 
and Mq — 800 GeV. Here the dotted line corresponds to the result computed without including the 
stop/stau coannihilation effects. 



So far, we have been focusing on the specific intermediate scale mirage mediation model 
with aj = Cj = 1 which might be obtained when all modular weights rii = 0. However as 
anticipated in the previous section, Oj = Cj = 1 is not necessarily a more favored choice 



than the other values of (aj,Ci) in Eq. (|3.2[) . Different values of and q, e.g. smaller but 
still non- negative values, are also equally plausible. Obviously, for a fixed value of Mq, the 
gaugino masses are not affected by changing Oj and q. However the low energy stop, stau 
and Higgs masses are somewhat sensitive to the values of Oj and Cj. They depend on Oj 
and Cj either through their boundary values at Mqut, or through their RG evolutions, or 
through the mass-mixing induced by the low energy A-parameters. 

The effects of changing ai and Cj on the RG evolution can be read off from the following 
one-loop RG equations for the Higgs and third generation sfermion mass-squares: 



ie 2 ^ 2 

dt ^" 

ic 2 ^ 2 

tc 2 2 

1 c 2 2 



3Xt-6gi\M2\' --gi\Mi\', 
5 

3Xb + Xr-6gl\M2f -^gi\Mi\^, 

5 



Xt + Xb 
2Xt 



^almi' - Gg^W - ^gl\Mi\^, 



fglm\'-%gl\Mi\\ 



167r2Am|3 = 2Xr -'^gl\Mi\\ (3.4) 

where 

Xt = 2yj{m\^ + + + Ajj^Q^^-J^ 

Xr = 2yl{ml^ + + ml.^ + ^1^,1.3^3). (3.5) 

These RG equations show that smaller Xi (I = t, r) increase the low energy soft mass- 
squares. Since aj and q determine the modulus- mediated trilinear A parameters and soft 
mass-squares at Mqut as Aij^ = (a^ + aj + ak)MQ and ihf = CiM^, smaller an = ch 
give smaller Xi without affecting the boundary values of squark and slepton masses at 
Mgut, eventually making the stop and stau masses at TeV scale larger. On the other 
hand, the consequence of smaller om = cm is more complicate as it depends on the relative 
importance of the Yukawa-induced RG evolution. It turns out that changing om = cm to 
smaller value makes the stop mass larger, while the stau mass smaller. 

In Figs, ^.b and Q.b, we depict the results for the case in which Oj = c, = 1/2 for 
both the matter and Higgs multiplets. As can be understood from the above discussion, 
this intermediate scale mirage mediation does not contain any parameter region of stop 
LSP, while having a larger parameter region of stau LSP (see Fig. |^.b). Another important 
feature is that the weak scale value of | m'jj^ \ becomes smaller compared to the case of 
Oi = Cj = 1, which is mainly due to smaller Xf. This results in smaller fi and rriA- 
Smaller fi makes the neutralino LSP have a sizable Higgsino component, while smaller ttia 
makes the pseudo-scalar resonance region disappear. Again the magenta region in Fig. 
|3[b corresponds to the parameter region giving the WMAP DM density (|3.3| ) under the 
assumption of pure thermal production. In this case, the neutralino pair annihilation into 
gauge boson pair becomes efficient due to the enhanced Higgsino component of neutralino 
LSP. Finally, the brown region is excluded by giving the Br(5 sj) smaller than the 
allowed range. 

Figs. |3|.c and Q.c are the result for the case with om = cm = 1 and an = ch = 0, 
while Figs, ^.d and ^.d are for the case with au = cm = 1/2 and an = ch = 0. The case 
of aM = Cm = 1/2 and an = ch = is quite similar to the case of Oj = Cj = 1/2: LSP 
is the lightest neutralino with a sizable Higgsino component for tan/3 < 20 (see Figs. |^.b 
and ^.d). On the other hand, the case of om = cm = 1 and oh = ch = is somewhat 
distinctive since there is no parameter region of stop or stau LSP and the WMAP DM 
density is obtained for a light neutralino mass m^o ~ 250 GeV, while in other cases the 
WMAP DM density is obtained for heavier m^o > 350 GeV. Again, the brown region is 
excluded by giving the Br (6 — > sj) smaller than the allowed range. 

3.2 Dark matter detections 

If neutralino LSP is the main component of the matter budget in the Milky Way, it might be 
detected through the elastic scattering with terrestrial nuclear target ^ . In the MSSM, 
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Figure 6: Scatter plot of the spin-independent neutralino-proton scattering cross section vs. 
in intermediate scale mirage mediation with {ai,Ci) specified in Eq. (3.2). 



t-channel Higgs boson and s-channel squark exchange processes contribute to the spin- 
independent (scalar) scattering between neutralino and nuclei. In many cases, dominant 
contribution to the scalar cross section comes from the Higgs exchange process which 
becomes bigger for larger tan/3, smaller Higgs masses, and mixed Bino-Higgsino LSP. 

In the specific intermediate scale mirage mediation model with = Cj = 1, the neu- 
tralino LSP is Bino-like and the mass of heavy CP-even Higgs boson is rather large when 
we require the neutralino to be LSP. It is thus expected that the elastic scattering cross 
section between neutralino DM and nuclei is rather small. In Fig. ^.a, we depict spin- 
independent (scalar) cross section asi of neutralino-proton scattering as a function of the 
LSP neutralino mass in this specific intermediate scale mirage mediation. Here we have im- 



posed the experimental bounds on the Higgs/sparticle masses and b ^ s'y branching ratio, 
and required that the hghtest neutrahno is the LSP. Red points in the figure correspond to 
the parameter values giving the WMAP DM density ( |3.3| ) under the assumption of pure 
thermal production, while the cyan points represent the parameter values for which the 
thermal production mechanism gives a smaller relic density. As expected, the cross section 
in the case of aj = Cj = 1 is quite small: asi ^ 5 x 10~^ pb, which is much smaller than 
the current experimental upper bound. It is even smaller than the sensitivity of future 
experiment such as SuperCDMS [^4| which would reach near 10~^ pb level. On the other 
hand, intermediate scale mirage mediations with different values of (aj,Cj) have a quite 
better prospect for direct detection. As can be seen from Fig. ^, most of the (red) WMAP 
points are above the sensitivity of SuperCDMS for the other three cases of different (aj, Cj). 
This is mainly due to the enhanced Higgsino component of the neutralino LSP and the 
reduced Higgs mass. 



Let us now examine gamma ray signals from DM annihilation in the galactic center, 
providing another feasible but indirect detection method for dark matter. The integrated 
gamma ray flux depends on the quantity J(Ar2), which is a measure of the cuspiness of 
the galactic halo density profile over a spherical region of solid angle AO,. In this paper, we 
use a conservative galactic halo model (isothermal halo density profile) which gives J ~ 30 
with the detector angular resolution AO = 10~^ sr and set -Ethr = 1 GeV for gamma 
ray energy threshold. Fig. |^ shows continuum gamma ray flux from the galactic center in 
intermediate scale mirage mediation scenarios under consideration, where red points give 
the WMAP value (|3.3| ) of the relic DM density. Here the four different choices of Oj = q do 
not lead to a dramatic difference in the gamma ray flux. The maximal value of flux given 
by the most favored WMAP (red) points is about few x 10^"'^^cm~^s^^ which is somewhat 
below the expected reach (~ 10~^'^cm~^s~^) of GLAST, although the (cyan) points giving 
smaller relic density can give a larger flux around 10~^''cm~^s~^. However, it should be 
noticed that our calculation for the gamma ray flux is based on a conservative halo density 
profile. If one uses an extreme halo model like the spiked profile the resulting gamma 
ray flux increases by a factor of ~ 10^. In this case, the gamma ray signals can be detected 
for a significant portion of the parameter space. A caveat is that the continuum gamma 
ray signals suffer from unknown astrophysical background. Recent observations of a bright 
gamma ray source in the direction of galaxy center by the Air Cherenkov Telescopes such 
as H.E.S.S. |3^] might be explained by an astrophysical process rather than the dark matter 
annihilation |37|. 



We finally notice an interesting enhancement of the gamma ray flux due to the Higgs 
resonance effect. Fig. ^ shows the gamma ray flux from the galactic center as a function 
of tan/3 in the specific intermediate scale mirage mediation with aj = Cj = 1 and Mq = 800 
GeV. One can see a clear enhancement of the flux around tan/3 ~ 22 for which rriA ~ 2m^. 
In this case, neutralino annihilation to heavy quarks is dominated and the subsequent 
quark hadronization produces many gamma rays. 
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Figure 7: Scatter plot of the continuum gamma ray flux vs. in intermediate scale mirage 
mediation with {ai,Ci) specified in Eq. (3.2). 



4. Neutralino DM for generic mirage messenger scale 

In the previous section, we have examined the prospect of neutralino DM in intermediate 
scale mirage mediation models (a = 1). As was discussed in section 2, in string compactifi- 
cations with non-trivial dilaton-modulus mixing, the anomaly to modulus mediation ratio 
a can have a more variety of values. In fact, the nature of neutralino LSP is somewhat 
sensitive to the value of a, typically it changes from Bino-like to Higgsino-like via Bino- 
Higgsino mixing region when a is increased from zero to a value of order unity. This feature 
is essentially due to the following behavior of the gaugino masses as a function of a: 



M3 : M2 : Ml ~ (1 - 0.3a)gl : (1 + 0.1a)g1 : (1 + 0.660)5?, 



(4.1) 
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Figure 8: Continuum gamma ray flux as a function of tan/? in the intermediate scale mirage 
mediation with Ui — Ci = 1 and Mq — 800 GeV. Note the resonant peak due to the psuedo-scalar 
Higgs resonance. 

If a increases from zero, the gluino mass decreases as M3 oc (1 — 0.3a). Smaller M3 then 
weakens the radiative electroweak symmetry breaking mechanism as it gives a smaller stop 
mass-square, thus leads to smaller |m//^ p and \fi\ at the weak scale. On the other hand, 
the Bino mass increases as Mi oc (1 + 0.66a), thus the lightest neutralino changes from 
Bino-like to Higgsino-like when a is varying from zero to a positive value of order unity. If 
a is further increased, eventually the model does not allow electroweak symmetry breaking. 
In this section, we extend the analysis of the previous section to the range of a from zero 
to the value at which the electroweak symmetry starts to be restored. 

4.1 Parameter region of neutralino LSP and thermal relic Density 

Again, let us first consider the case with aj = q = 1. We will treat Mg and a as free 
parameters, while focusing on tan/3 = 10 and 35. Figs, ^.a and ^.b show how some of the 
superparticle masses vary as a function of a for a fixed Mq = 800 GeV. For a < 1, the LSP 
is the lightest neutralino which is mostly Bino, and thus its mass varies as m^o ~ Mi oc 
(1+0. 66a). In the range of 1 < a < 1.8, stau or stop becomes the LSP. For 1.8 < a < 2, the 
lightest neutralino which is now mostly Higgsino becomes the LSP. If a increases further, 
the model does not allow electroweak symmetry breaking. 

In Figs. |9|.c and ^.d, the two distinct magenta regions seperated by stop/stau LSP 
regions give the WMAP DM density, 0.085 < QoMh'^ < 0.119, under the assumption 
that all neutralino DMs are produced by the conventional thermal production mechanism. 
Below (above) these magenta regions, Q^h'^ < 0.085 (> 0.119). In the Bino-like LSP 
region, stop-neutralino coannihialtion plays a crucial role to get the WMAP DM density for 
tan/3 = 10, while stau-neutralino coannihilation or pseudoscalar Higgs resonance processes 
are important for tan/3 = 35. For Higgsino-like LSP, the charged Higgsino xf ^-iid two 
neutral Higgsinos Xi)X2 nearly degenerate. Then the dominant annihilation processes 
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Figure 9: Sparticle masses vs. a for tan/3 = 10 and tan/3 = 35 in case with = c, = 1. The lower 
figures show the parameter space of neutralino LSP and its thermal relic density on the plane of 
(a, Mo). 



are the neutralino pair annihilation into gauge bosons, and the neutralino-chargino co- 
annihilation into fermion pair ||3^. These annihilations of Higgsino-like LSP are very 
efficient, so that the relic mass density is too small unless is quite heavy. Indeed, from 
Figs. |9|.c and ^.d, we can see that the WMAP DM density is obtained only for Mq > 2.2 
TeV in the Higgsino LSP region around a ~ 1.8. However it should be stressed that the 
cyan regions of Figs, ^.c and |9|.d can be allowed if some part of DM were produced by 
non-thermal mechanism. Such parameter region contains a ~ 1.8 and Mq ~ 1 TeV for 
which the neutral Higgsino with m^o ~ 200 GeV is the LSP and the stop is rather light as 
~ 250 GeV. The brown regions are excluded by the 6 ^ 57 constraint. On the brown 
region in Fig. 9.c, the chargino loop contribution to 6 — > 57 dominates, which results in 
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Br(6 — > ,37) smaller than the experimentally allowed range. 





Let us now consider the case with Oj = c, = 1/2. Obviously, smaller (a,, Cj) give smaller 
stop/stau mass-squares at Mqut- However, as was anticipated in the previous section, Xi 
(I = t, b, r) which govern the RG evolution of stop/stau mass-squares (see Eq. |3.4| ) become 
smaller also, which would increase the stop/stau masses at the weak scale. Together with 
the reduction of Ah^q^Us^ ^^^^ effect on the RG evolution eventually makes the physical 
lighter stop mass m^^ larger compared to the case with aj = q = 1. On the other hand, 
stau masses are more affected by the change of the boundary values, thus their weak scale 
values become lighter compared to the case of = Cj = 1. Smaller Xt leads to also a 
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Figure 11: The results for the case in which an = ch = and om = cm = 1- 



smaller at the weak scale, resulting the reduction of the Higgsino mass fj, and the 

pseudoscalar Higgs boson mass m^. Figs. 10. a and 10. b show all of these features. Again, 
as a increases, the neutralino LSP changes from Bino-like to Higgsino-like. Comparing to 
Fig. 1^, the lighter stop becomes heavier, while the lighter stau and the pseudoscalar Higgs 
become lighter. As a consequence, the stop LSP region disappears, but the stau LSP region 
at large tan/? becomes larger. The magenta regions of Figs, ^.c and 10. d correspond to 
the parameter region giving the WMAP DM density under the assumption of pure thermal 
production. They clearly show the Higgsino-like LSP at a > 1 and also the pseudoscalar 
Higgs resonance effect for the Bino-like LSP at smaller a. 



-23- 





Figure 12: The results for the case in which an — cu — and cm = cm = 1/2. 



Fig. 11 shows the results for the case in which oh = ch = for the Higgs multiplets, 
while Cm = cm = 1 for the quark/lepton matter multiplets. A characteristic feature of this 
case is that the lightest neutrahno is the LSP over the entire region of parameter space 
allowing the electroweak symmetry breaking. Compared to the case in which Oj = q = 1 
for both the Higgs and matter multiplets, Xj {I = t, b, r) for the RG evolution ( |3.4| ) have 
smaller values, while the boundary values of stop/stau mass-squares remain the same. This 
results in heavier stop and stau at the weak scale. Except for the absence of stop/stau 
LSP region, other features are somewhat similar to other cases. The brown regions are 
excluded by the b ^ sj constraint. On the brown region with small Mq in Fig. 11. c, 
the chargino loop contribution to 6 — > 57 dominates, which results in Br(6 — > 57) smaller 
than the experimentally allowed range. On the other hand, the charged Higgs boson loop 
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becomes significant in the large Mq region, making the predicted Br(6 57) exceed the 
experimental bound. The region between those two brown regions is allowed due to the 
cancellation between the chargino and charged Higgs boson loop contributions. 



Finally, Fig. 12 is for the case with an = ch = and an = cm = 1/2. The results are 
quite similar to the case in which Oj = Cj = 1/2 for both the Higgs and matter multiplets. 
The a = 2 region of this case corresponds to the TeV scale mirage mediation model 
proposed in as a model to minimize the fine tuning for the electroweak symmetry 
breaking in the MSSM. 

4.2 Dark matter detections 
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Figure 13: Spin-independent neutralino and proton scattering cross section in case with 
1. 



To see the prospect of direct DM detection, spin-independent cross section of the 
neutralino-proton scattering is presented in Fig. ^ for the case with Oj = Cj = 1. Here, we 
imposed the experimental bounds on sparticle and Higgs masses, and 6 — > 57 branching 
ratio. In the figures, the red points give the WMAP DM density: 0.085 < fl^h'^ < 0.119, 
the cyan corresponds to the region giving O.^.h'^ < 0.085, and the rest gives ^x^'^ ^ 0.119, 
under the assumption of pure thermal production of neutralino DM. One can notice that 
there are two distinct branches of the WMAP points which correspond to the Bino branch 
and the Higgsino branch, respectively. In our scan, Higgsino-like LSP gives a larger asi 
for a given m^- The dominant contribution to asi usually comes from the Higgs exchange 
process which becomes significant if the LSP neutralino is a mixed Bino-Higgsino state. On 
the other hand, for aj = q = 1, the LSP neutralino is either Bino-like or Higgsino-like since 
the mixed Bino-Higgsino region gives a stop or stau LSP. Therefore, it is expected that the 
cross section is rather small for the case with Oj = Cj = 1. Indeed, Fig. 13 shows that the 



predicted values are all less than the current and near future experimental sensitivity. 
However, the prospect of direct DM detection is dramatically changed if one considers 



other choices of Oj and Cj. Fig. 14 shows the predictions for spin-independent cross section of 
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Figure 14: Spin-independent neutralino and 
(ai, Ci) giving a mixed Bino-Higgsino LSP over 



proton scattering cross section for other values of 
a significant fraction of the parameter space. 



the neutralino-proton scattering for the three other choices of (a,, Cj) giving a mixed Bino- 
Higgsino LSP over a significant fraction of the parameter space and also a reduced value 
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of the pseudoscalar Higgs mass. These values of a, and Cj give heavier stop, thereby the 
6 ^ 57 constraint becomes less significant compared to the case with = Cj = 1. Again, 
the red points represent the parameter values giving the WMAP DM density 0.085 < 
n^h^ < 0.119, the cyan points give ^i^/i^ < 0.085, and the rest stands for ^i^h^ > 0.119, 
under the assumption of thermal production of neutralino LSP. As expected, the scattering 
cross sections are largely enhanced compared to the case with Oj = q = 1. Now, much 
of the WMAP points give asi exceeding the sensitivity limit of the planned SuperCDMS 
experiment. If one includes the cyan points, the cross section can be much bigger, reaching 
even at the current CDMS sensitivity limit. 
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Figure 15: Continuum (a and b) and monochromatic (c and d) gamma ray flux from the Galactic 
Center vs. for the case with ai = Ci = 1 and tan/3 = 10 or 35. 



Gamma rays induced by neutralino annihilation in Galactic Center might provide an 



indirect detection of neutralino DM. Fig. 15 shows the predicted continuum (a and b) 
and monochromatic (c and d) gamma ray fluxes from the Galactic Center as a function 
of the LSP neutralino mass for the case with Oj = c, = 1. Here we chose the same halo 
density profile as the previous section, giving J(AJ7 = 10~^sr) ~ 30. The red points in 



the figures give the WMAP DM density, while the cyan and the rest give Q^h'^ < 0.085 
and ^x^'^ ^ 0.119, respectively, under the assumption of pure thermal production. Again 
the WMAP points have two distinct branches, the Bino-branch and the Higgsino-branch. 
Including the cyan points giving smaller thermal relic DM density, the case with Oj = q can 
give a continuum gamma ray flux up to 2 x lO^^^cm^^s^^ and lO^^'^cm^^s^^ for tan /? = 10 
and 35, respectively. This maximum flux of the continuum gamma rays barely touch the 
expected reach of GLAST. However, the real gamma ray flux can be much bigger than 
these predictions if the actual halo density profile is denser than the assumed profile. For 
Higgsino LSP, unsuppressed annihilation into W or Z boson pair is the major source of 
continuum gamma rays. As can be noticed from Fig. |l5|, for some parameter values, the 
gamma ray flux from Bino LSP is largely enhanced by the pseudoscalar Higgs resonance 
effect. 
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Figure 16: Continuum (a and b) and monochromatic (c and d) gamma ray flux from the Galactic 
Center vs. for the case in which an — cn — 0, om = cm — 1/2, and tan/3 — 10 or 35. 



In fact, it is quite nontrivial to discriminate the continuum gamma rays produced by 
neutralino annihilation from the diffuse galactic gamma ray backgrounds. On the other 



hand, the monochromatic gamma ray hne from XX ~^ 77 or 1^ can be considered as a 
'smoking gun' signal of WIMP dark matter. Figs. |l^.c and 15. d show the gamma ray 
hne flux produced by neutralino pair annihilation in Galactic Center for the case with 
Oj = Cj. One can notice that there is a clear distinction between the Bino and Higgsino 
LSP regions. The gamma ray line flux ranges from 10^^^ to 10~^^cm~^s~^ for the Bino 
LSP branch of WMAP points, while it ranges from 10" 1^ to 10" 1^ cm s for the Higgsino 
LSP branch. For Higgsino-like LSP, the gamma ray line flux comes dominantly from the 
W^Xi loop diagrams resulting in a large cross section for xX ~^ 11 or 7Z |3^. While 
GLAST will probe the photon energies only up to 300 GeV with a low energy threshold. 
Atmospheric Cherenkov Telescopes (ACT) such as H.E.S.S. will be able to cover higher 
photon energy ranges and probe the gamma ray flux down to 10"^^cm~^s"^. The predicted 
monochromatic fluxes in Figs. |l^.c and ^.d are still below this sensitivity limit. However, 
as we have stressed, these results are based on a rather conservative halo density profile. In 
view of that the predicted flux can increase even by a factor of 10^ if one uses an extreme 
halo model like the spiked profile, the monochromatic gamma ray signal for the Higgsino 
dark matter might be measurable in case of a cuspy halo density profile. 

As we have anticipated, other values of (ai,Ci) specified in ( |3.2D allow a mixed Bino- 
Higgsino LSP over a significant fraction of parameter space. It is thus expected that those 
other cases can give a larger gamma ray flux compared to the case with Oj = Cj = 1. In 
Fig. |l6|, we depicted the results for the case with an = ch = and om = cm = 1/2. 
Indeed, this case gives a larger fiux, although not dramatically different. The red WMAP 
points can give a continuum gamma ray flux up to 3 x 10"^^cm~^s~^, while the cyan points 
giving smaller thermal relic DM density can reach up to 2 x 10"^'^cm"^s"^. The maximal 
flux of monochromatic gamma ray is about 10"^^cm~^s"^ for the red WMAP points and 
about 7 X 10~^^cm"^s~^ for the cyan points. Again, these results are obtained for the 
conservative halo density proflle giving J(Af] = 10~^sr) ~ 30. The real gamma ray flux 
can be significantly bigger than these predictions if the actual halo density profile is denser 
than the assumed proflle. 



5. Conclusions 



In this paper, we have examined the prospect of neutralino dark matter in mirage media- 
tion scenario of SUSY breaking in which soft masses receive comparable contributions from 
modulus mediation and anomaly mediation. Depending upon the model parameters, espe- 
cially the anomaly to modulus mediation ratio, the nature of the lightest neutralino changes 
from Bino-like to Higgsino-like via Bino-Higgsino mixing region. For Bino-like LSP, the 
conventional thermal production mechanism can give a right amount of relic DM density, 
i.e. the WMAP observation 0.085 < QoAih'^ < 0.119, through the stop/stau- neutralino 
coannihilation process or the pseudo-scalar Higgs resonance effect. In overall, compared 
to the mSUGRA scenario, a signiflcantly larger fraction of the parameter space can give 
the WMAP DM density under the assumption of thermal production, while satisfying all 
known phenomenological constraints. This is partly because the lightest neutralino is a 
mixed Bino-Higgsino over a sizable fraction of the parameter space. 



We also studied the detection possibilities of neutralino dark matter in mirage medi- 
ation. For the parameter region giving the WMAP density of Bino-like or Higgsino-like 
LSP, direct detection via elastic scattering between neutralino DM and nuclear target turns 
out to be mostly under the sensitivity of near future experiments. However the other pa- 
rameter region giving the WMAP density of mixed Bino-Higgsino LSP predicts typically 
a cross section above the expected sensitivity limit of SuperCDMS. The continuum and 
monochromatic gamma ray fluxes from neutralino annihilation in Galactic Center have 
been analyzed also. Generically, Higgsino-like LSP gives a larger gamma ray flux than 
Bino-like LSP, however the continuum gamma ray flux from Bino LSP can be significantly 
enhanced for some particular parameter values due to the pseudo-scalar Higgs resonance 
effect. Although the gamma ray fluxes predicted within a conservative halo model are 
below the sensitivity of ongoing and planned experiments, it might be detectable if the 
actual halo density is denser than the conservative profile used in our analysis. 
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Appendix A. 

In this appendix, we summarize the notations and conventions used in this paper. The 
quantum effective action in N = 1 superspace is given by 



/ 
/ 




) 




(5.1) 



where the gauge kinetic terms are written as a Z?-term operator to accommodate the 
radiative corrections to gauge coupUngs, and the elhpsis stands for the irrelevant higher 
dimensional operators. The Kahler potential K is expanded as 

K = Ko{Ta, TX) + Zi{TA, ri)$*e2^'"^''$i + (5.2) 

where Va and <l>j denote the visible gauge and matter superfields given by 

= + V2 9i>' + e^F\ 



= -Oa'^eA''^ - io'^ex" + ie^ex"" + -e'^e^n", (5.3) 



1 
2 

and Ta = (C, T) are the SUSY breaking messengers including the conformal compensator 
superfield C = Co + O'^F^ and the modulus superfield T = Tq + y/2ef + The 
radiative corrections due to renormalizable gauge and Yukawa interactions can be encoded 
in the matter Kahler metric Zi and the gauge coupling superfield Ga which is given by 

Ga = Re(/a) + AG,, (5.4) 

where fa is the holomorphic gauge kinetic function and AGa includes the T^-dependent 
radiative correction to gauge coupling. The superpotential is expanded as 

W = Wo{T) + ^Xijk{T)<^i^j<^k + (5.5) 

where Wo{T) is the modulus superpotential stabilizing T. Here we do not specify the mech- 
anism to generate the MSSM Higgs parameters fi and B, and treat them as free parameters 
constrained only by the electroweak symmetry breaking condition. For a discussion of ^ 
and B in mirage mediation, see Ref. [^. 

For the canonically normalized component fields, the above superspace action gives the 
following form of the running gauge and Yukawa couplings, the super symmetric gaugino- 
matter fermion coupling Cx^p, and the soft SUSY breaking terms: 

^ = Re{Ga), Vijk = -^=^===, 
9a Ve '^oZiZjZk 

Cx^ = iV2 {cl)*T''AX'' - X'^T^ckii,,) , 

>Csoft = -mU'<p'* - ^M.A'^A'' + ^Ai,fcy,,fc<^>J0'= + h.c. j , (5.6) 



where 



Ma = F^aAln(Re(G„)), 



for 



^e-K^Z,ZjZk, 

mj = -F^F^*dAdB In (e-^o/^^i) (5.7) 
F^ = -e^°/^{dTdT*)-\DTWoy, 

F^ = m*y^ + \dTKoF'^ ("^3/2 = e^"/'Wo). (5.8) 



In the approximation ignoring the ofF-diagonal components of Wij = YlpqVipgyjpq^ the 1- 
loop RG evolution of soft parameters is determined by 



167r 



167r' 



dln/x 

2 dAijk 
din /J, 

2 dm^ 



= 2 



3tr(T2(Adj))+^tr(T2(</)^)) 

i 



+ 



+ 



k 



3,k 



(5.9) 



where the quadratic Casimir C2{4>i) = {N"^ — 1)/2N for a fundamental representation 
of the gauge group SU{N), C2{(t>i) = qf for the U{1) charge Qi of ^j. 

In mirage mediation, soft terms at Mqut are determined by the modulus mediation 
of 0{F'^ /T) and the anomaly mediation of 0{F'-'' /Stt'^Cq) which are comparable to each 
other. In the presence of the axionic shift symmetry 

U{1)t ■■ Im(T) + real constant (5.10) 

which is broken by the non-perturbative term in the modulus superpotential 

Wo = w- ^e-"^, (5.11) 

one can always make that TO3/2 and are simultaneously real. Also since F^ /T ~ 
m3/2/47r2, we have 



(5.12) 



Then, upon ignoring the parts of 0{F'^ /Stt'^T), the resulting soft parameters at McuT are 
given by 



Ma = Mo + — ^ 6. -2 



167r2 
'3/2 



m,- = m 



~2 "^3/2 , . „ /"'<V2\2 . 



(5.13) 



where 



Mo = F^aTlnRe(/„), 
Aijk = (a^ + a^- + afc)Mo = F'^ drHer^'^ ZiZjZk), 
c^M^ = -\F'^\^dTdf\vL{e-^'''^Zi), 



rhf = 



(5.14) 



o 
o 

DO 

o 
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o 
o 
o 
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and 



ba = -3tr (T„2(Adj)) + J^tr (Tf (<^,)) , 

i 

7. = 25^5^C7f(</,,)-^^b,,feP, 



jk 



ii = Svr 



a jk 

2 dji 



(5.15) 



where coij = '^kiUikiU^kl assumed to be diagonal. Here we have used that \jk are 
T-independent constant as ensured by the axionic shift symmetry U{1)t- 

Let us now summarize our conventions for the MSSM. The superpotential of canoni- 
cally normaUzed matter superfields is given by 



W = yoHd ■ QD^ + VLHd ■ LE^ - yuH^ ■ QU" - fiH^ ■ Hu, 



(5.16) 



where the SU{2)l product is H ■ Q = eabH"'Q^ with ei2 = —€21 = 1, and color indices are 
suppressed. Then the chargino and neutralino mass matrices are given by 



where 



Mc 



M 



N 



-M2 92{H^)] 

/ -Ml -j=,9Ym j=,9Ym \ 

-M2 j^92{H'd) -^92{H^) 



(5.17) 



(5.18) 



in the field basis 



Vi'^^ = -i [b, w\ iHl mi) , 



) 



(5.19) 



for = {W^^iW^)/V2. 

The one-loop beta function coefficients ba and anomalous dimension 7j in the MSSM 
are given by 



63 = -3, 62 = 1, 

iHu = ^92 + \9y - 



o 
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DO 

o 
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o 
o 
o 
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(5.20) 

where g2 and = ^/3/5gl denote the SU{2)l and U{1)y gauge couplings. The 9i and ji 
which determine the soft scalar masses at Mqut are given by 

^Hu = 3^1 +9y- ^VtiaHr. + (iQz + au-,), 

dHd = ^92 +9y - ^vli^'Hi + «Q3 + - "^yliO'Hi + flLa + ^^3) 

^Qa = y^i + 3^1 + ^9y - '^{vtiO'Hu + aQs + aua) + vliaHa + 0Q3 + or>3))<^3a, 

^ 16 2 16 2 . 2/ N c- 

^t/a = + -^9y - %t + «Q3 + ac/3)'^3a, 

/I 16 2 4 2 , 2/ Nf 
^Da = y53 + g5y - 42/^ (a^d + aq^ + aDgj^sa, 

^La = 3^2 +9y - '^vliO'Hi + aL3 + aE3)ha, 

^Ea = - %r(aj^d + «i3 + a£;3)'^3a, (5.21) 



and 



3 4 11 



7i/„ = 2^2 + Y^y - 3% 6y„ 



where 



3 4 11 4 o 27 27 

= 2^2 + yfy - 3^6 - y^V' 
3 11 

iQa = -8^3 + 25'2 + Y^S'y - {vIK + ybK)^3a, 
88 

7[/a = -853 + y5y - '^'VtbytSsa, 
22 

7D„ = -85I + —gy - 2ylby^S3a, 

3 4 11 4 27 c 

7L„ = 2f2 + Y^y - Vrbyrha, 

jE^ = 22g^ - 2yXj3a, (5.22) 



V = -y53 - 3^2 - Y^y + 6% + Vb, 

K = ~ Y^^ ~ ~ ^* 

byr = -332' - 35?^ + 3yl + Ayl (5.23) 
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